Dec. 2015 航 天 控制 


Vol. 33, No.6 Aerospace Control “4 


DOI:10.16804/j.cnki.issn1006-3242.2015.06.007 


太阳 同步 轨道 立方 星 任务 轨 妃 演 化 分 析 


雷 博 持 ” 李 明 涛 : 
1. 中 国 空间 技术 研究 院 西 安 分 院 , 西 安 710100 
2. 中 国 科 学 院 空间 科学 与 应 用 研究 中 心 ,北京 100190 


摘 要 立方 星 是 近年 来 微小 卫星 领域 的 研究 热点 , 文章 以 太阳 同步 轨道 
( SSO) 立方 星 任务 为 例 ,研究 地 球 非 球形 引力 项 J, 和 大 气 阻 力 对 立方 星 轨 道 演 
化 的 影响 。 分 别 研究 了 轨道 高 度 、 太 阳 活 动 指数 F10.7 和 大 气 模型 对 SSO 立方 
星 轨道 演化 的 影响 ,进行 了 数值 仿真 验证 。 研 究 表 明 : 大 气 阻 力 摄 动 对 SSO 立方 
星 轨 道 寿命 影响 很 大 ,影响 结果 与 轨道 高 度 之 间 存 在 非 线 性 关系 ; 太阳 活动 高 、 
低 年 时 的 轨道 演化 差异 明显 , 且 演 化 量 随 F10.7 的 变化 是 非 线性 的 ; 不 同 大 气 模 
型 下 的 轨道 演化 情况 存在 差异 ,大 多 数 模型 所 对 应 的 轨道 演化 量 的 量 级 基本 一 致 。 
上 述 影响 分 析 对 SSO 立方 星 任务 轨道 的 选 定 及 轨道 控制 具有 一 定 的 参考 价值 。 
关键 词 ”立方 星 ; 轨道 演化 ; 大 气 模 型 ; 太阳 活动 指数 
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人 Abstract Cube-Sat is the research focus in the field of micro satellite in recent years. In this paper, mis— 
sion that with a cubessat in the Sun synchronous orbit ( SSO) is taken as an example, then the orbit evolu— 
tion of SSO cubesat is studied considering both the influence of atmosphere drag and that of J, coefficient 
( the Farth’ s oblateness perturbations to the second-order zonal coefficient). Orbit altitude, solar activity 
index and atmosphere models are researched respectively to show their impacts on the evolution of orbit, and 
numerical simulations are carried out. The work above indicates some conclusions as follows: perturbation 
force that produced by air drag plays an important role to decide the lifetime of SSO cubessat, and there is a 
nonlinear relationship between the orbit altitude and the magnitude of orbital decay which is induced by atmos— 
phere drag. Orbital decay based on high solar activity is obviously different from that based on low solar activity, 
and the relationship between F10.7 and the magnitude of orbital change is not linear. Considering various air 
drag models separately, the orbit evolutions are not equal , but their order of magnitude are almost the same. 
The analysis above is helpful to both the orbital choosing of SSO CubeSat and its orbit control. 
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立方 星 这 一 概念 最 早 由 美国 加 州 理 工大 学 和 斯 
坦 福 大 学 在 1999 年 提出 ,立方 星 是 标准 化 的 微小 
卫星 ,尺寸 10 em x 10 cm x 10 cm 质量 约 1 kg ,为 一 
个 基本 单元 1 U。 以 此 为 基础 可 以 有 2 U、3 U 等 立方 
星 。 立 方 星 具 有 很 多 特点 ,不 仅 成 本 低 , 开发 时 间 也 
较 短 ,易于 实现 标准 化 、 模 块 化 ,易于 技术 更 新 "-” 。 

近年 来 ,利用 多 卫星 系统 的 空间 任务 越 来 越 多 ， 
借助 标准 化 的 微小 卫星 来 实现 分 布 式 空间 系统 极 具 
发 展 前 景 。 分 布 式 空间 任务 的 出 现 , 拓展 了 立方 星 
在 科学 领域 的 应 用 。 目 前 已 有 越 来 越 多 的 立方 星 被 
送 上 太空 。 国 外 的 立方 星 应 用 主要 分 为 3 个 方面 : 
以 空间 科学 和 教学 为 核心 目标 ,如 荷兰 代 尔 夫 特 理 
工大 学 的 Delfi-c3; 以 载荷 任务 为 目标 ,如 QuakeSat 
搭载 了 地 震 电磁 探测 载荷 ; 也 有 以 新 技术 验证 为 目 
标的 立方 星 任务 ,如 加 拿 大 的 先进 纳 星 太空 实验 室 
CANX 卫星 系列 ”。 国 内 很 多 高 校 及 科研 机 构 也 在 
积极 开展 立方 星 的 相关 研究 。 

体积 小 、 重 量 轻 的 自身 特点 决定 了 立方 星 大 多 
没有 在 轨 机 动能 力 , 考虑 无 主动 控制 时 立方 星 轨道 
的 自然 演化 就 成 为 一 个 重要 研究 方向 。 太 阳 同 步 轨 
道 ( SSO) 具有 星 上 光照 条 件 好 、 过 同一 纬度 的 当地 
时 间 相 同等 特点 ,已 被 广泛 应 用 于 各 种 空间 任务 。 
因此 ,本 文选 取 了 SSO 立方 星 任务 进行 轨道 演化 研 
究 , 考 虑 地 球 非 球形 引力 项 J, 和 大 气 阻力 摄 动 的 影 
响 , 以 平均 半 长 轴 随 时 间 的 变化 作为 摄 动 效应 的 指 
标 。 以 轨道 分 析 软 件 为 仿真 平台 ,先进 行 了 摄 动力 
简化 以 及 立方 星 面 质 比 对 轨道 演化 的 影响 分 析 , 以 
确立 数值 仿真 的 前 提 条 件 ; 接着 研究 了 轨道 高 度 、 太 
阳 活 动 指 数 和 大 气 模型 这 3 个 方面 对 SSO 立方 星 
轨道 演化 的 影响 。 最 终 得 出 结论 。 


1 卫星 轨道 动力 学 模型 


1.1 坐标 系 定 义 

研究 地 球 非 球形 摄 动 及 大 气 阻 力 摄 动 对 轨道 演 
化 的 影响 时 ,涉及 到 地 心 惯性 坐标 系 和 地 心地 固 坐 
标 系 ,如 图 1 所 示 。Oxsyiz ( 标记 为 Si) 为 地 心 惯 性 
坐标 系 ，Ox。y,z。( 标记 为 Se) 为 地 心地 固 坐标 系 ，a 
为 Greenwich 赤 经 。 

地 心 惯性 坐标 系 Si : 坐标 系 原点 在 地 心 , Z 轴 
垂直 于 地 球 赤道 平面 ; XY 轴 在 地 球 赤道 平面 内 , X 
轴 指 向 春分 点 位 置 ; 了 轴 满 足 右 手 螺旋 定 则 。 动 力 
学 方程 一 般 建立 在 惯性 坐标 系 中 。 


1 地 心 惯 性 坐标 系 和 地 心地 固 坐 标 系 


地 心地 固 坐标 系 Se : 坐标 系 原 点 在 地 心 , Z 轴 
垂直 于 地 球 赤道 平面 ; XY 轴 在 地 球 赤道 平面 内 , X 
轴 指 向 Greenwich 子午 圈 与 赤道 交点 ; 了 轴 满 足 右 
手 螺旋 定 则 。 地 球 非 球形 摄 动 一 般 在 地 心地 固 坐标 
系 中 表示 。 

1.2 动力 学 模型 

在 惯性 坐标 系 下 , 卫星 动力 学 模型 可 用 式 ( 1) 

表示 。 


7 =- Kr+f (1) 


式 中 : 7 为 地 心 距 矢量 , w 是 地 心 引力 常数 ,7 为 地 
心 距 大 小 , f 为 地 球 中 心 引力 之 外 的 所 有 摄 动力 的 
合力 和 包括 地 球 非 球 形 引 力 、 大 气 阻力 、 太 阳光 压 、 
第 三 体 引 力 摄 动 等 。 

将 式 (1) 在 惯性 坐标 系 下 展开 并 进行 积分 , 可 
以 得 到 任意 时 刻 卫 星 在 惯性 坐标 系 下 的 运动 状态 。 


2 地球 非 球形 摄 动 和 大 气 阻力 摄 动 分 析 


对 于 低 轨 道 来 说 ,卫星 所 受到 的 摄 动力 主要 是 
地 球 非 球形 引力 摄 动 、 大 气 阻力 摄 动 、 太 阳光 压 摄 动 
以 及 第 三 体 引 力 摄 动 。 设 有 表示 地 球 中 心 引力 加 
速度 的 大 小 ,上 表示 其 他 类 型 的 摄 动 加 速度 的 大 小 。 
与 地 球 中 心 引力 相 比 ,用 % = 大力 来 表征 各 摄 动力 
的 量 级 估算 值 ,如 表 1 所 示 。 其 中 大 气 阻力 摄 动 与 
第 三 体 引 力 摄 动 力 均 为 10” ,但 第 三 体 引 力 摄 动 要 
小 于 大 和 气 阻 力 摄 动 , 且 轨 道 高 度 越 低 , 大 气 阻力 摄 动 
的 作用 越 大 。 由 此 可 见 ,对 于 低 轨 卫星 , 六 摄 动 最 
重要 ,其 次 是 大 气 阻力 振动 ”。 我 们 可 以 对 摄 动力 
进行 简化 , 仪 考虑 ., 摄 动 和 大 气 阻力 摄 动 的 影响 。 
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表 1 各 种 摄 动 力 的 量 级 
摄 动力 类 型 ”地球 非 球形 ” 太阳 光 压 大 气 阻力 ”第 三 体 引 力 
量 级 ( 7 ) 10-3 102? 1077 1077 


2.1 /项 摄 动 

在 两 体 问 题 中 ,认为 地 球 是 一 个 均匀 球体 ,卫星 
运行 在 地 球 中 心 引力 场 中 。 但 实际 上 地 球形 状 是 不 
规则 性 的 ,质量 分 布 也 不 均匀 , 常 将 地 球 引 力 势 范 数 
表示 成 如 下 形式 : 


CR 
VrsgA) = 全 + 全 忆 ( 放 ViPsing+ 全 ， 


ES 
> ( 一 ) "( Cm CoSmA + Sim sinmA) Psing (2) 

n=2 m=1 三 
式 中 : 7 为 到 观测 点 的 地 心 距 , 由 为 地 心经 度 , A 为 


地 心 纬 度 , Rs 是 赤道 半径 ，P,，P,, 是 正则 勒 让 德 
多 项 式 , C,, ,5,, ,J 为 引力 场 系数 。 

对 引力 势 范 数 (2) 求 梯度 , 可 以 得 到 卫星 在 地 
心地 固 坐标 系 下 的 加 速度 , 再 通过 坐标 转换 转 到 地 
心 惯性 坐标 系 , 即 可 以 得 到 惯性 系 下 卫星 的 动力 学 
方程 。J, 项 的 影响 要 比 其 他 项 大 100 倍 以 上 , 通常 
仅 考 虑 J 项 摄 动 的 影响 。 

2.2 大 气 阻力 摄 动 

大 气 阻力 是 非 保守 力 , 在 其 影响 下 ,卫星 轨道 会 
发 生 衰减 。 大 气 阻力 摄 动 对 卫星 所 产生 的 加 速度 如 
下 式 : 

fy = -二 Co 二 pf， 区 
V, =v -Vv, (3) 
式 中 : f 为 大 气 阻 力 摄 动 加 速度 , V, 为 卫星 相对 大 
气 的 相对 速度 矢量 , 到 为 卫星 相对 大 气 的 相对 速度 
的 大 小 ,v 和 v, 分 别 为 卫星 和 大 气相 对 地 心 惯性 坐 
标 系 的 速度 矢量 , p 为 大 气 密度 , C, 为 阻力 系数 。 
S/m 为 面 质 比 参数 。 

另外 , 分析 大 气 阻 力 摄 动 时 , 必须 考虑 以 下 几 方 
面 因素 : 

1) Cn ;阻力 系数 一 般 与 卫星 表面 材料 有 关 , 常 
取 2.2; 

2) S/m , 面 质 比 与 有 效 迎 风 面 积 有 关 , 即 要 考 
虑 卫星 姿态 和 形状 ,通常 选用 一 个 等 效 面 质 比 来 进 
行 计算 ; 

3) 7?。, 大气 转速 比较 复杂 , 约 为 地 球 自转 速度 
的 0.8 ~1.4 倍 ,通常 认为 两 者 相等 ; 

4) p ;大 气 密度 是 极其 复杂 的 问题 ,与 轨道 高 
度 、 温 度 、 太 阳 活 动 指数 等 密切 相关 , 变化 幅 值 非常 


大 。 已 有 的 大 气 模型 也 很 多 , 如 指数 模型 , Jacchia 
77 模型 、NRLMSISE2000 模型 等 , 其 中 NRLM- 
SISE2000 大 气 模型 是 由 美国 海军 研究 实验 室 于 
2000 年 在 MSISE-90 模型 的 基础 上 发 展 而 出 ,是 最 
新 更 新 的 一 个 大 气 模 型 ,与 其 他 模型 相 比 具有 很 多 
优势 。 
2.3 ”项 摄 动 和 大 气 摄 动 影响 下 的 卫星 运动 方程 
基于 2.1 节 和 2.2 节 的 分 析 , 可 以 得 到 同时 考 
虑 /项 摄 动 和 大 气 阻力 摄 动 时 ,卫星 在 惯性 系 下 的 
运动 方程 如 下 式 : 


3 
A ee 
r r T r 
1 5S 
SC 
y My., 3,Rs,? 
7 =- 和 所- 和 插 J 
r T r 
1 5S 
SC Dh 
3 E\2 z 
i ) (全 
r FE 
1 5 
Fe (4) 


给 定 卫星 初始 时 刻 的 运动 状态 后 , 对 式 ( 4) 进 
行 数值 积分 ,就 可 以 得 到 在 J, 摄 动 和 大 和 气 阻 力 摄 动 
同时 影响 下 ,卫星 轨道 的 演化 进程 ”。 


3 仿真 分 析 


在 轨道 仿真 软件 STK 中 建立 一 个 SSO 轨道 ,将 
其 偏心 率 、 近 地 点 幅 角 、 真 近 点 角 设 置 为 0, 降 交点 
地 方 时 设置 为 18 时。 改变 轨道 高 度 、 设 定 轨道 预报 
器 就 可 以 得 到 满足 文章 要 求 的 各 种 仿真 条 件 。 

首先 分 析 了 面 质 比 对 轨道 演化 的 影响 并 进 
行 仿真 ,以 确立 出 基本 的 仿真 初始 条 件 。 然 后 分 
别 仿真 分 析 了 轨道 高 度 、 太 阳 活 动 指数 以 及 不 同 
大 气 模型 对 轨道 演化 的 影响 。 由 于 瞬时 轨道 根 
数 波动 性 比较 大 ,无 法 表现 卫星 受 摄 运动 的 本 
质 ,文中 用 平均 轨道 根 数 下 的 半 长 轴 作 为 摄 动 效 
应 的 衡量 指标 。 
3.1 面 质 比 的 选 定 

由 2.2 节 可 知 ,在 研究 大 气 阻力 摄 动 时 , 必须 考 
虑 卫星 面 质 比 的 影响 ,本 文 研究 对 象 是 1 U 立方 星 ， 
边 长 为 10 cm, 质量 约 为 1 kg, 所 以 面 质 比 的 变化 范 
围 为 :0.01 ~0.014。 在 STK 中 ,考虑 .项 和 大 气 模 
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型 为 NRLMSISE2000 时 的 大 气 摄 动 , 在 0.01 ~ 6969 
0.014 之 间 以 0.0005 为 步 长 改变 面 质 比 ,仿真 不 同 ER 
面 质 比 对 500km SSO 立方 星 轨道 演化 的 影响 ,利用 Re 
所 得 9 个 数据 点 绘图 ,如 图 2 所 示 。 横 坐标 是 面 质 
比 , 纵 坐标 代表 半 长 轴 的 衰减 量 。 
24 针 6966 
23 69065.5 
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2 SS0 立方 星 轨 道 演化 与 面 质 比 的 关系 


从 图 2 可 以 看 出 , 面 质 比 对 轨道 演化 有 很 大 影 
响 , 仿真 分 析 时 必须 合理 选择 面 质 比 ; 轨道 演化 情况 
大 致 与 面 质 比 成 线性 关系 。 文 中 以 下 章节 的 仿真 
中 , 面 质 比 默认 选 定 为 0.01m /kg。 
3.2 ”轨道 高 度 对 轨道 演化 的 影响 

以 NRLMSISE2000 大 气 模 型 为 例 ,研究 大 气 对 
不 同 高 度 SSO 立方 星 轨道 演化 的 影响 ,仿真 时 间 
180d。NRLMSISE2000 大 气 模型 中 , 太阳 活动 指数 
F10.7 按 默认 值 为 130。 在 300 ~ 700km 范围 内 以 
100km 为 步 长 改变 轨道 高 度 , 所 得 到 的 半 长 轴 演 化 
情况 汇总 如 表 2 所 示 , 从 仿真 结果 中 选 出 400km， 
600km 高 度 轨道 的 演化 图 ,分 别 如 图 3 和 4 所 示 。 
由 此 可 以 得 到 180d 内 轨道 衰减 的 情况 。 图 3 中 ,在 
大 气 阻 力作 用 下 ,145d 左右 时 轨道 高 度 直 线 下 降 ， 
卫星 陨落 ,寿命 终止 。 图 4 中 , 仿真 结束 时 ,轨道 衰 
减 了 约 3.6km。 
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表 2 轨道 高 度 对 轨道 演化 的 影响 ( F10.7 =150) 


轨道 高 度 /km 300 400 500 600 700 
半年 内 半 长 硅 命 等 命 
16 3.6 1 
轴 变 化 /km 23d 145d 


从 表 2 可 以 看 出 :1) 大 气 阻 力 摄 动 对 轨道 演化 
( 半 长 轴 衰 减 ) 影响 很 大 , 且 轨 道 高 度 越 低 影响 越 明 
显 ; 2) 高 度 300km 的 轨道 , 寿命 约 23d4。 高 度 为 
400km 时 ,轨道 寿命 约 150d。 而 对 于 500km 高 的 轨 
道 , 半 长 轴 衰 减 了 16km。 由 此 可 见 , 轨 道 的 演化 情 
况 与 轨道 高 度 间 存在 非 线 性 关系 。 
3.3 ”太阳 活动 对 轨道 演化 的 影响 

太阳 活动 指数 F10. 7 指 太 阳 发 出 的 波长 为 
10.7cm 的 电磁 辐射 强度 , 其 变化 没有 严格 规律 ,一 
般 由 观测 得 到 。 

F10.7 的 长 期 预测 具有 不 准确 性 ,以 2015 年 6 
月 1 日 为 例 ,STK 轨道 分 析 软 件 8.1.1 根据 2007 年 
5 月 24 日 的 观测 数据 ,预测 的 F10.7 参数 为 77. 7; 
STK 轨道 分 析 软 件 9.2.1 根据 2009 年 4 月 29 日 的 
观测 数据 ,预测 的 F10.7 参数 为 101.4; NOAA 根据 
2012 年 11 月 26 日 的 观测 数据 ,预测 的 F10.7 参数 
为 117.2”。 由 此 可 见 ,实现 以 年 为 单位 的 F10. 7 
参数 预报 , 误差 很 大 。 应 该 采用 尽 可 能 新 的 预测 
数据 。 

由 于 了 10.7 的 不 确定 性 , 导致 大 气 密度 相应 地 
也 存在 不 确定 性 。 因 此 ,在 考虑 大 气 阻 力 摄 动 时 ,有 
必要 研究 太阳 活动 指数 F10.7 对 轨道 演化 的 影响 。 

根据 太阳 活动 指数 随 年 份 变化 的 趋势 ,选择 太 
阳 活 动 高 、 低 年 时 的 F10. 7 参数 分 别 为 200 和 
70 。 然 后 分 别 分 析 太 阳 活 动 高 、 低 年 时 , 大气 对 
不 同 高 度 SSO 演化 情况 的 影响 。 
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当 轨 道 高 度 为 300 ~ 1000km 时 ,可 以 分 别 得 到 
太阳 活动 高 、 低 年 时 轨道 半 长 轴 的 演化 情况 ,汇总 如 
表 3。 


表 3 F10.7 参数 对 轨道 演化 的 影响 


轨道 高 “太阳 活动 高 年 时 半 长 轴 太阳 活动 低 年 时 半 长 轴 
度 /km ”变化 (F10.7 =200) 变化 ( F10.7 =70) 

300 寿命 16d 寿命 50d 

400 寿命 95d 18 km 

500 37 km 2 km 

600 9 km 300 m 

700 2.5 km 140 m 

800 1 km 基本 保持 不 变 

900 300 m / 

1000 基本 保持 不 变 / 


选取 太阳 活动 低 、 高 年 时 500km SSO 的 演化 情 
况 ,分 别 如 图 5 和 6 所 示 。 由 图 5 和 6 可 以 看 出 ,对 
于 500km 的 SSO 轨道 ,太阳 活动 低 、 高 年 时 的 轨道 
演化 情况 差距 很 大 ,仿真 结束 时 ,轨道 半 长 轴 分 别 豪 
减 了 2km 和 37km。 
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图 5 太阳 活动 低 年 时 的 轨道 演化 情况 
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图 6 太阳 活动 高 年 时 的 轨道 演化 


从 表 3 可 以 看 出 , 对 于 同样 高 度 的 SSO 立方 
星 , 太 阳 活 动 高 、 低 年 时 的 轨道 演化 差别 很 大 ; 太阳 
活动 指数 越 高 ,可 以 忽略 大 气 阻 力 摄 动 的 最 低 轨道 
高 度 也 越 高 。 一 般 认 为 180 天 内 半 长 轴 衰 减 量 小 于 
1 km 时 ,就 可 以 忽略 大 气 阻力 的 影响 。 

下 面 以 500 km SSO 轨道 为 例 , 研究 轨道 高 度 不 
变 时 , 半 长 轴 衰 减 量 随 F10.7 的 变化 ,以 得 到 F10.7 
对 轨道 衰减 的 影响 曲线 图 。 在 70 ~ 200 之 间 改 变 
F10.7, 步 长 取 10。 利 用 所 得 数据 点 绘图 ,如 图 7 所 
示 。 由 图 7 可 知 ,对 同一 高 度 的 轨道 ,太阳 活动 指数 
越 高 ,轨道 衰减 就 越 快 ; 轨道 衰减 量 与 F10.7 呈现 出 
非 线性 关系 。 
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7 半 长 轴 衰减 量 与 F10.7 之 间 的 关系 
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3.4 ”大 气 模型 对 轨道 演化 的 影响 

由 2.2 节 可 知 , 分析 大 气 阻力 摄 动 时 , 必须 考虑 
大 气 模 型 的 影响 。 目 前 已 经 有 多 个 大 气 模型 ,然而 
这 些 模型 的 统计 精度 只 有 15% 左右 。 

以 高 度 为 500 km 的 典型 SSO 轨道 为 例 , 选用 
不 同 的 大 气 模型 ,对 轨道 半 长 轴 随 时 间 的 演化 进行 
分 析 , 仿真 初始 条 件 与 上 文 相同 。 分 别 选用 了 8 个 
大 气 模 型 : Cira 72 模型 , Exponential 下 arth 模型 , Har- 
ris-Priester 模型 , Jacchia-Roberts 模型 ,Jacchia-1971 
模型 ,MSIS 1986 模型 ,US Standard Atmosphere 模型 
以 及 NRLMSISE 2000 模型 。 各 模型 中 涉及 到 太阳 
活动 指数 F10.7 的 ,都 按照 默认 设置 为 150。 

对 于 指数 大 气 模 型 ,仿真 时 要 输入 参考 高 度 和 
标准 高 度 , 按 式 (5) 计算 ”。 


(rn) (5) 


式 中 : 万 为 标准 高 度 , 有 为 参考 高 度 , + 为 卫星 质心 
的 地 心 距 。 = 0.1 ( 常 取 jw < 0.2)，, HH =37.4 km， 
ro。 = 有 +6371 km。 因 此 ,对 于 500 km 轨道 ,计算 可 
得 参考 高 度 是 37.4 km ,标准 高 度 是 60. 887 km。 


H=H,+ 
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分 别 设置 不 同 的 轨道 预报 器 ,以 改变 大 气 模型 。 
仿真 180d, 得 到 不 同 大 气 模型 下 半 长 轴 的 演化 情 
况 , 如 表 4 所 示 。 


表 4 不 同 大 气 模型 对 轨道 演化 的 影响 ( 180d) 


大 气 模型 半 长 轴 变 化 /km 
Cira 72 18 
Exponential— Earth 3 
Harris Priester 反光 
Jachia Roberts 18 
Jacchi 1971 18 
MSIS 1986 16 
NRLMSISE2000 16 
US Standard 10 


从 表 4 可 以 看 出 ,在 不 同 大 气 模型 的 影响 下 , 轨 
道 演化 情况 不 尽 相 同 ,大 多 数 模型 的 影响 效果 的 量 
级 约 为 16 km; 与 其 他 模型 相 比 , 指数 大 气 模型 影响 
下 的 半 长 轴 演 化 情况 悬殊 ,反映 出 指数 大 气 模型 的 
精度 不 够 高 。 


4 结论 


立方 星 是 近年 来 纳 卫星 领域 的 研究 热点 。 本 文 
对 SSO 立方 星 任务 的 轨道 演化 情况 进行 了 仿真 分 
析 , 得 出 如 下 结论 : 

1) 对 于 300 ~700 km 轨道 高 度 的 立方 体 卫 星 ， 
大 气 阻 力 摄 动 对 轨道 寿命 有 明显 影响 ,导致 半 长 轴 
不 断 衰减 ,轨道 越 低 影 响 越 大 , 且 影 响 结 果 与 轨道 寿 
命 成 非 线 性 关系 ; 

2) 太阳 活动 对 大 气 阻 力 摄 动 有 显著 影响 , 太阳 
活动 高 、 低 年 时 的 轨道 演化 情况 差距 很 大 。 在 分 析 
大 气 阻力 摄 动 时 ,必须 考虑 F10.7 参数 的 影响 , 且 影 
响 效果 与 F10.7 间 存 在 非 线性 关系 ; 

3) 选取 不 同 大 气 模型 时 ,轨道 演化 情况 不 完全 
相同 。 大 多 数 大 气 模型 之 间 的 轨道 演化 情况 差距 并 
不 大 ,指数 模型 与 其 他 模型 相 比 存在 一 定 偏差 。 

综 上 所 述 ,考虑 轨道 高 度 、 太 阳 活 动 及 大 气 模型 
的 耦合 影响 , 对 单 颗 SSO 立方 星 任务 轨道 的 选 定 具 
有 参考 意义 。 另 外 , 单 星 的 轨道 演化 是 研究 立方 星 
编队 的 基础 , 星 间 无 碰撞 的 安全 飞行 高 度 的 确定 、 星 


间 相 对 距离 的 演化 情况 及 编队 构 形 控制 方法 可 以 作 
为 下 一 步 的 研究 方向 。 
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